Hydrodynamic limit for weakly asymmetric 
simple exclusion processes in crystal lattices 

Ryokichi Tanaka 
Abstract 

We investigate tiie hydrodynamic limit for weakly asymmetric simple exclu- 
sion processes in crystal lattices. We construct a suitable scaling limit by using a 
discrete harmonic map. As we shall observe, the quasi-linear parabolic equation 
in the limit is defined on a flat torus and depends on both the local structure of 
the crystal lattice and the discrete harmonic map. We formulate the local ergodic 
theorem on the crystal lattice by introducing the notion of local function bundle, 
which is a family of local functions on the configuration space. The ideas and 
methods are taken from the discrete geometric analysis to these problems. Results 
we obtain are extensions of ones by Kipnis, 011a and Varadhan to crystal lattices. 

1 Introduction 

The purpose of this paper is to discuss the hydrodynamic limit for interacting particle 
systems in the crystal lattice. Problems of the hydrodynamic limit have been studied 
intensively in the case where the underlying space is the Euclidean lattice. We extend 
problems to the case where the underlying space has geometric structures: the crystal 
lattice. The crystal lattice is a generalization of classical lattice, the square lattice, the 
triangular lattice, the hexagonal lattice, the Kagome lattice (Figure[T]i and the diamond 
lattice. Before explaining difficulties for this extension and entering into details, we 
motivate to study these problems. 

There are many problems on the scaling limit of interacting particle systems, which 
have their origins in the statistical mechanics and the hydrodynamics. (See HJ, Iil5l 
and references therein.) The hydrodynamic Hmit for the exclusion process is one of 
the most studied models in this context. Here we give only one example for exclusion 
processes in the integer lattice, which is a prototype of our results, due to Kipnis, Olla 
and Varadhan (H]). From the view point of physics and mathematics, it is natural to ask 
for the scaling limit of interacting particle systems evolving in more general spaces and 
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Figure 1 : Crystal Lattices 

to discuss the relationship between macroscopic behaviors of particles and geometric 
structures of the underlying spaces. In this paper, we deal with the crystal lattice, which 
is the simplest extension of the Euclidean lattice Z"'. Although the crystal lattice has 
periodic global structures, it has inhomogeneous local structures. 

On the other hand, crystal lattices have been studied in view of discrete geometric 
analysis by Kotani and Sunada (Q, IfTOll . IfTTII . and the expository article |fT6l ). They 
formulate a crystal lattice as an abelian covering graph, and then they study random 
walks on crystal lattices and discuss the relationship between asymptotic behaviors of 
random walks and geometric structures of crystal lattices. In ifTOl . they introduce the 
standard realization, which is a discrete harmonic map from a crystal lattice into a 
Euclidean space to characterize an equilibrium configuration of crystals. In |9|, they 
discuss the relationship between the Albanese metric which is introduced into the Eu- 
clidean space, associated with the standard realization and the central limit theorem 
for random walks on the crystal lattice. Considering exclusion processes on the crystal 
lattice, one is interested to ask what geometric structures appear in the case where the 
interactions depend on the local structures. 



2 



Given a graph, the exclusion process on it describes the following dynamics: Parti- 
cles attempt to jump to nearest neighbor sites, however, they are forbidden to jump to 
sites which other particles have already occupied. So, particles are able to jump to near- 
est neighbor vacant sites. Then, the problem of the hydrodynamic limit is to capture 
the collective behavior of particles via the scaling limit. If we take a suitable scaling 
limit of space and time, then we observe that the density of particles is governed by a 
partial differential equation as a macroscopic model. Here it is necessary to construct 
a suitable scaling limit for a graph and to know some analytic properties of the limit 
space. 

A crystal lattice is defined as an infinite graph X which admits a free action of a 
free abelian group F with a finite quotient graph Xq. We construct a scaling limit of 
a crystal lattice as follows: Let be a positive integer Take a finite index subgroup 
A^F of F, which is isomorphic to NZ'' when F is isomorphic to Z'^. Then we take the 
quotient of X by A^F-action: X^. We call this finite quotient graph Xi^i the N -scaling 
finite graph. The quotient group '■- TjNT acts freely on Xf^. Here we consider 
exclusion processes on X^. To observe these processes in the continuous space, we 
embed X^ into a torus. We construct an embedding map Oat from X^ into a torus by 
using a harmonic map O in the discrete sense in order that the image (^M(Xf^) converges 
to a torus as goes to the infinity. (Here the convergence of metric spaces is verified 
by using the Gromov-Hausdorff topology, however, we do not need this notion in this 
paper) Then we obtain exclusion processes embedded by Oat into the torus. 

In this paper, we deal with the simplest case among exclusion processes: the sym- 
metric simple exclusion process and its perturbation: the weakly asymmetric simple 
exclusion process. In the latter case, we obtain a heat equation with nonlinear drift 
terms on torus as the limit of process of empirical density (Theorerr l3.1l and Examples 
below). We observe that the diffusion coefficient matrices and nonlinear drift terms can 
be computed by data of a finite quotient graph and a harmonic map O. (See also 
examples in Section l272l ) The hydrodynamic limit for these processes on the crystal 
lattice is obtained as an extension of the one on Z''. So, first, we review the outline 
of the proof for Z"', following the method by Guo, Papanicolaou and Varadhan in ||5|. 
Since the lattice Z'' is naturally embedded into R'', the combinatorial Laplacian on the 
scaled discrete torus converges to the Laplacian on the torus according to this natural 
embedding. The local ergodic theorem is the key step of the proof since it enables us to 
replace local averages by global averages and to verify the derivation of the limit par- 
tial differential equation. It is formulated by using local functions on the configuration 
space and the shift action on the discrete torus. The proof of the local ergodic theorem 
is based on the one-block estimate and the two-blocks estimate. Roughly speaking, the 
one-block estimate is interpreted as the local law of large numbers and the two-blocks 
estimate is interpreted as the asymptotic independence of two different local laws of 
large numbers. 

Second, we look over the outline of the proof for the crystal lattice. There are two 
main points with regard to the difference between Z'' and the crystal lattice, that are 
the convergence of the Laplacian and the local ergodic theorem. Although the crystal 
lattice X is embedded into an Euclidean space by a harmonic map <1), the combinatorial 
Laplacian on the image of the A^-scaling finite graph <^i^{Xii) does not converge to the 
Laplacian on the torus straightforwardly. It is proved by averaging each fundamental 
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domain by T-action because of the local inhomogeneity of the crystal lattice. Thus, it 
is necessary to obtain the local ergodic theorem compatible with the convergence of 
the Laplacian. Furthermore, it is also necessary to obtain the local ergodic theorem 
compatible with the local inhomogeneity of the crystal lattice. For these reasons, we 
have to modify the local ergodic theorem in the case of crystal lattices. To formulate 
the local ergodic theorem in the crystal lattice, we introduce the notion of T-periodic 
local function bundles. A F-periodic local function bundle is a family of local functions 
on the configuration space which is parametrized by vertices periodically. Moreover, 
we introduce two different ways to take local averages of a F-periodic local function 
bundle. The first one is to take averages per each fundamental domain as a unit. The 
second one is to take averages on each F-orbit. The local ergodic theorem in the crystal 
lattice is formulated by using F-periodic local function bundles, two types of local av- 
erages and the F^y-action on the A^-scaling finite graph X^- In fact, we use only special 
F-periodic local function bundles to handle the weakly asymmetric simple exclusion 
process. The proof of this local ergodic theorem is also based on the one-block esti- 
mate and the two-blocks estimate. Proofs of these two estimates are analogous to the 
case of the discrete torus since we use the fact that the whole crystal lattice is covered 
by the F-action of a fundamental domain in the first type of the local average and we 
restrict to a F-orbit in the second type of the local average. In this paper, we call the 
local ergodic theorem the replacement theorem and prove it in the form of the super 
exponential estimate. The derivation of the hydrodynamic equation is the same manner 
as the case of the discrete torus. 

Let us mention related works. Interacting particle systems are categorized into the 
gradient system and the non-gradient system, according to types of interactions. We 
call the system the gradient system when the interaction term is represented by the 
difference of local functions. Otherwise, we call the system the non-gradient system. 
We mention a recent work on the non-gradient system by Sasada lfT3l . The symmetric 
simple exclusion process is a model of the gradient system. Our problems essentially 
correspond to problems for the gradient system since the hydrodynamic limit for the 
weakly asymmetric simple exclusion process is reduced to the one for the symmetric 
simple exclusion process, following 1 8 1 . As for the hydrodynamic limit on spaces other 
than the Euclidean lattice, Jara investigates the hydrodynamic limit for zero-range pro- 
cesses in the Sierpinski gasket (16 1). As for the crystal lattice, there is another type of 
the scaling limit. In 1 14|, Shubin and Sunada study lattice vibrations of crystal lattices 
and calculate one of the thermodynamic quantities: the specific heat. They derive the 
equation of motion by taking the continuum limit of the crystal lattice. As a further 
problem, we mention the following problem: Recently, attentions have been payed for 
interacting particle systems evolving in random environments (e.g., [JJ, [3J and lilSl ). 
For example, the quenched invariance principle for the random walk on the infinite 
cluster of supercritical percolation of Z'' with > 2 is proved by Berger and Biskup 
(UJ). Their argument is based on a harmonic embedding of percolation cluster into 
R''. Our use of the harmonic map O and local function bundles will play a role in the 
systematic treatment of particle systems in more general random graphs. Furthermore, 
the hydrodynamic limit on the inhomogeneous crystal lattice is considered as the case 
where the crystal lattice has topological defects. This problem would be interesting in 
connection with material sciences. 
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This paper is organized as follows: In Section |2] we introduce the crystal lattice 
and construct the scaling limit by using discrete harmonic maps. In Section[3] we for- 
mulate the weakly asymmetric simple exclusion process on the crystal lattice and state 
the main theorem (Theorem |TT). In Section|4] we introduce F-periodic local function 
bundles and show the replacement theorem (Theorem |4.H . We prove the one-block es- 
timate and the two-blocks estimate. In Section |5] we derive the quasi-linear parabolic 
equation, applying the replacement theorem and complete the proof of Theorem 13.11 
Section |6] is Appendix; A. We prove some lemmas related to approximation by com- 
binatorial metrics to complete the scaling limit argument. Section |7] is Appendix;B. 
We refer an energy estimate of a weak solution and a uniqueness result for the partial 
differential equation to this appendix. 

Landau asymptotic notation. Throughout the paper, we use the notation a - o^ Xo 
mean that a — > as — » oo. We also use the notation a = to mean that a — > as 



2 The crystal lattice and the harmonic realization 

In this section, we introduce the crystal lattice as an infinite graph and its realization 
into the Euclidean space. 

2.1 Crystal lattices 

Let X - (V, E) be a locally finite connected graph, where V is a set of vertices and E 
a set of all oriented edges. The graph X may have loops and multiple edges. For an 
oriented edge e e £, we denote by oe the origin of e, by te the terminus and by e the 
inverse edge of e. Here we regard X as a weighted graph, whose weight functions on 
V and E are all equal to one. 

We call a locally finite connected graph X = (V, E) a T-crystal lattice if a free 
abelian group F acts freely on X and the quotient graph r\X is a finite graph Xq - 
(Vo, £■())■ More precisely, each cr e F defines a graph isomorphism cr : X X and the 
graph isomorphism is fixed point-free except for cr - id. In other words, a F-crystal lat- 
tice X is an abelian covering graph of a finite graph Xq whose covering transformation 
group is F. 

2.2 Harmonic maps 

Let us construct an embedding of a F-crystal lattice X into the Euclidean space R"' of 
dimension d - rankT. 

Given an injective homomorphism ^ : F — > R'' such that there exits a basis 
Ml, . . . , Mrf e R'', 




then we define a harmonic map associated with (p. 
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Definition 2.1. Fix an injective homomorphism (p as above. We call an embedding 
C) : X — > R"', a (p-periodic harmonic map if^ satisfies the fallowings: O is Y-periodic, 
i.e., far any x & V and any cr e F, 0(cr;ic) — 0(x) + cp{cr) and O is harmonic, i.e., far 
any x eV, Iiee£j<l*(fe) - <l)(oe)] = 0, where E_, - {e e E \ oe - x). 

We note that a 0-periodic harmonic map O depends on and call O a periodic 
harmonic map in short when we fix some (p. 

Fore e £0, we take a lift e e £ of e, and define v(e) := <i>{te)-<i>(oe) e R"'. By the F- 
periodicity, v(e) does not depend on the choices of lifts. For v(e) = (v'l (e), . . . , v,/(e)) e 
W', let us define adx t/-matrix by 



Here the matrix is symmetric and positive definite. We call the matrix D the diffusion 
coefficient matrix. 
Examples 

0. The one dimensional standard lattice. 

Oa. The one dimensional standard lattice X which we identify the set of vertices V 
with Z and the set of (unoriented) edges with the set of pairs of vertices {{x,x + 
1) I X e Z). Now Z acts freely on X by the additive operation in Z and the quotient 
finite graph consists of one vertex and one loop as un oriented graph. When we 
regard X as an oriented graph, we add both oriented edges to X and the quotient 
graph consists of one vertex and two oriented loops (Figure |2]i. 

Let us choose a canonical injective homomorphism : Z — > R. In our formula- 
tion, choose a basis ei = 1 in R and define : Z — > R by setting </>(«) = ne\ for 
n e Z so that 0(Z) = {nei | n e Z). By identifying the set of vertices of X with Z, 
we define an embedding map = (p{x), x e Z. This embedding map O is a 
Z-periodic harmonic map. In this case, D = 1/2. 

Ob. Let us give another example of periodic harmonic map for the one dimensional 
standard lattice X. Now we define a Z-action on X in the following way: For 
cr e Z, X e y, define crx :- Icr + x. Then this induces a free Z-action on X 
and the quotient graph consists of two vertices and two edges between them as 
an unoriented graph. Let (p be the injective homomorphism as the same as in 
Example Oa. We define an embedding map O : X — > R by setting 0(0-0) :- 
-1- 0(cr), 0(0-1) -\ + 4>(cr). Then O is a periodic harmonic map. The image 
of O is not isomorphic to the previous one (Figure|3). In this case, D = 5/4. 

1. The square lattice. 

1 a. The square lattice has the standard embedding in R^ and this embedding is shown 
to be periodic and harmonic in our sense in the following. We identify the set of 
vertices of the square lattice X with 7? and the set of edges with the set of pairs 
of vertices {(x, x + (1,0)), (x, x + (0, 1)) | x e I?)]. Now Z^ acts freely on X by 
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Figure 2: The one dimensional standard lattice and the quotient graph in Example Oa. 



Figure 3: The image of <I) in Example Ob. 

the additive operation in 7} and the quotient graph is the bouquet graph with one 
vertex and two unoriented loops. When we regard X as an oriented graph, we 
add both oriented edges to X and the quotient finite graph is the bouquet graph 
with one vertex and four oriented loops. 

Let us choose a canonical injective homomorphism : — > R^. That is, 
choose a basis \e\ - (1,0), - (0, 1)) in and define : ^ R^ by setting 
^{{m,n)) - me\ + ne^ for (m, n) G l? so that 0(Z^) - \Y!i=\ I integers). 
By identifying the set of vertices of X with 7?, we define an embedding map 
0(x) = X € 7} . This embedding map O is a Z^ -periodic harmonic map. In 
this case, D = (1/2 ^O^j ^ 

lb. Let us give another example of periodic harmonic map for the square lattice X. 
Choose a basis {mi = (1,0), = (1 /2, 1)) in R^ and define : Z^ -> R^y set- 
ting 0((m, n)) - mu\ + nu^ for (m, ri) G l? so that 0(Z^) - {2?=i ^/M; I integers). 
In the same way as above Example la, we define an embedding map <l>(x) - <p(x), 
X G Z^. (Figure |4]) This embedding map O is a Z^ -periodic harmonic map. In 
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Figure 4: The square lattice embedded as in Example lb. and the quotient graph. 

Ic. Let us give an example an embedding map <i> which is periodic but not harmonic. 
We choose an action of on the square lattice X in the following way: Again, 
we identify the set of vertices V of X with Z^. For cr = (0-1,0-2) € Z^, x - 
(xi,X2) e V, define crx :- (cri + x\,2cr2 + X2). Then this induces a free Z- 
action and the quotient graph consists of two vertices, two edges between them 
and one loop on each vertex (two loops) as an unoriented graph. Let be the 
same as in Example la. We define an embedding map O' : X — > by setting 
O'(o-(0,0)) = (0,0) + 0(cr), a)'(cr(0, 1)) = (1, 1/2) + (^(o-) for o- e 7?. Then O' 
is periodic but not harmonic since for x - (0, 0) e 7}, YjeeE I'^'it^) ~ •5'(oe)] = 
(1, 0) + (-1, 0) + (1, 1/2) + (1, -1/2) = (2, 0) 9i (0, 0). 

2. The hexagonal lattice. 

The hexagonal lattice admits a free Z^-action with the quotient graph consist- 
ing of two vertices and three edges as an unoriented graph. We define a fun- 
damental subgraph D by setting the set of vertices {xo,x\,X2,x-i] and the set 
of (unoriented) edges {{xq, x\), (xq, X2), (xq, X3)). Then the hexagonal lattice has 
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a subgraph isomorphic to D and is covered by copies of the subgraph trans- 
lated by the Z^-action. Choose a basis [ui = ( Vs, 0), U2 = ( V3/2, 3/2)) in 
and define : ^ R" by setting (f){(m,n)) = mui + nu2 for (m,n) G Z^ 
so that (f>(Z^) - ^iUi I ki integers). We define an embedding map O by 

setting <l)(crxo) = (0,0) + 0(cr), (D(crxi) = (-V3/2, 1/2) + 0(cr), (t(o-X2) = 
( V3/2, 1/2) + (^(o-) and <^{o-x^) = (0, -1) + 0(cr) for cr e 1?. (Figureg]) Then O 

is a periodic harmonic map. In this case, D = 



3/8 
3/8 





Figure 5: The hexagonal lattice embedded as in Example 2. and the quotient graph. 

3. The Kagome lattice. 

The Kagome lattice admits a free Z^-action with the quotient graph consisting of 
three vertices and six edges (two edges between each pair of vertices) as an un- 
oriented graph. We define a fundamental subgraph D by setting the set of vertices 
{xo, x\,X2, ^3, xn] and the set of (unoriented) edges {(xq, x\), (xq, X2), (xq, X3), (xq, X4)). 
Then the Kagome lattice has a subgraph isomorphic to D and is covered by 
copies of the subgraph translated by the Z^-action. Choose a basis {mi = ( V3, 0), M2 — 
(V3/2,3/2)) in R^ and define : Z^ ^ R^ as the same as in Example 2. 
We define an embedding map Q> by setting <l)(crxo) - (0,0) + <;*(<t), <l)(crxi) = 
(- V3/2,O)+0(cr), (D(o-X2) = (- V3/4,-3/4)+</.(o-), (D(o-X3) = ( V3/2, 0)+0(cr), 
0(crx4) = ( V^/4, 3/4) + 0(cr) for cr € 7?. (Figure ID) Then (D is a periodic har- 
monic map. In this case, D - 

Remark. 



3/8 
3/8 
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Figure 6: The Kagome lattice embedded as in Example 3. and the quotient finite graph. 

The notion of periodic harmonic map on F-crystal lattice is studied by Kotani and 
Sunada and including the standard realization which they introduced in 1 10] as a special 
case. They use harmonic maps to characterize equilibrium configurations of crystals. 
In fact, a periodic harmonic map is characterized by a critical map for some discrete 
analogue of energy functional. The standard realization is not only a critical map but 
also the map whose energy itself is minimized by changing flat metrics on torus with 
fixed volume. (More precisely, see|10|). The existence of periodic harmonic map 
for every injective homomorphism producing lattices in R'' and the uniqueness up to 
translation is proved in Theorem 2.3 and Theorem 2.4 in llTOl . 

2.3 Scaling Limits 

Let us construct the scaling limit of the crystal lattice. Suppose that T is isomorphic 
to Z^. Let > 1 be an arbitrary positive integer and A^F the subgroup isomorphic to 
NZ'^. The subgroup A^F acts also freely on X and its quotient graph Nr\X is also a 
finite graph Xn = (Vn, Em). Then Y/NF s Z^'/NZ'' acts freely on X^. We call Xn the 
N -scaling finite graph. The map 

— O : X ^ W', 
N 
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satisfies that (l/N)0(o-^x) - (l/N)(b(x) + ij/(cr) for all x e V and all cr e T since O is 
F-equivariant. We have the torus T*^ := equipped with the flat metric induced 

from the Euclidean metric. The torus depends on iff, however, we do not specify it in 
the following. Then the map (1/A^)<1> : X — » R'' induces the map 

We call the N-scaling map. (Figur^) 




N 



Figure 7: The image of the A^-scaling finite graph by a harmonic map in the covering 
space 

Next, we observe convergence of the combinatorial Laplacian on Xff. Since the 
degrees of x e Vn might be different, depending on each x, we consider "average" of 
the combinatorial Laplacian on a fundamental domain. 

Let d(x) be the degree of a vertex x e Vn, i.e., the cardinality of the set E^^j^ = {e e 
En \ oe - x]. Define the combinatorial Laplacian associated with - {Vn,En) 
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acting on the space of continuous functions C('F') by 

A^7((Da,(x)) := y [mN(te)) - 7(<Div(oe))J , 

^^^^ eeE^.. 

for J € C('r') and x € Vjv. We show that the combinatorial Laplacian converges to the 
Laplacian on T'' in the following sense: For every twice continuous derivative functions 
/ G C^(T''), for every x G T^, for each x g Vq, take an arbitrary sequence of vertices 
xn G such that x^v is a Uft of x and 0;v(-^iv) — > x as A'^ — > oo, then by the Taylor 
formula, 

( d 



^ ' V 1=1 ' 1J=1 ' 



N{xN))vi{e)vj{e) 



■ On- 



Since O is harmonic, 
1 



, xeVaeeEfi^,^ i,j=l J 



Since YjxsVo Ile€£„,;c„ Vi(e)v^(e) = HeeSo v;(e)v/e), the last term is equal to 2VDVy(x), 
where D is a diffusion coefficient matrix and VDV = and D = 

{dij)i,j=\,...4- 



3 Hydrodynamic limit for exclusion processes 

We formulate the symmetric simple exclusion process and the weakly asymmetric sim- 
ple exclusion process in crystal lattices. As we see below, the former is a particular case 
of the latter. 

Let Xm = (Vat, E^) be the A^-scaling finite graph of X. We denote the configuration 
space by Ziv := {0, 1)^'". We denote the configuration space for the whole crystal lattice 
X - (y, £) by Z := {0, 1)^. Each configuration is defined by t; = {rix)xeVN ^ with 
//x = or 1 and by 77 G Z in the same way. 

We consider the Bernoulli measure and Vp on Z^v, Z, respectively, for < p < 1 . 
They are defined as the product measures of the Bemoulli measure vj, on {0, 1), where 
yi(0)= l-p,yi(l)=p. 

Let L^{Zn, v^) be the -space of R-valued functions on Z^. The action of Tf^ on 
Xf^ hfts on Z/v by setting {cnf)x := T]a-^x for o" g F^r and x G Vff. The group F^r also acts 
on L\Zn,v^) by crFiT]) := Fia-^rj) for F g L\Zn,v^). For e g Ejy and 77 g Zn, we 
denote by 77^ the configuration defined by exchanging the values of rjoe and ?7,e. i.e.. 



nx--- 



T]te X = Oe 
Voe X — tC 

rjx otherwise. 
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For each e € E^, we define the operator : L^(Zff,v^) L?(Zff,v^) by setting 
n^Firj) :- F{jf ) - F(ri). We see that 77"^ = rj'' and n^F - n-^F for e € Ef^. The above 
notations also indicate corresponding ones for Z = {0, 1)*' the configuration space on 
the whole crystal lattice. 

The symmetric simple exclusion process is defined by the generator acting on 
LHZm,v^) as 

LNF{r]) = ^ 2 n^Fir]), F e L\Zn, 

The weakly asymmetric simple exclusion process is defined as a perturbation of 
the symmetric simple exclusion process. We denote by C''^([0, T] x T') the space of 
continuous functions with continuous derivatives in [0, T] and the twice continuous 
derivatives in T"'. For each function H € C''^([0, T] x T''), the weakly asymmetric 
simple exclusion process on is defined by the generator acting on L^(Zn,v^) as 



eeE, 



where 



c"(e, T], t) rioe (1 - 71,^} exp [H(t, ^N{te)) - H(t, Oiv(oe))] . 

Here Oat : — > T'' is the A^-scaling map. The meaning of the perturbation is as 
follows: We introduce a "small" drift depending on space and time in particles. In the 
original process, a particle jumps with rate 1 /2 from oe to te {e is an edge) at time t, 
while in the perturbed process, a particle jumps approximately with rate 

Therefore, the external field which is now (1 /2N)'VH gives a small asymmetry of the 
order 1 /N in the jump rate. Notice that we obtain the symmetric simple exclusion 
process when H is constant. 
Remark 

The weakly asymmetric simple exclusion process which we introduced here does 
not include the well-studied case where for one dimensional lattice, the external field is 
(l/2N)E for some constant E > Q and its limit equation produces the viscous Burgers 
equation (e.g., |l2l). This process corresponds to the case with VH = E which we do 
not treat here. 

Let Z)([0, T~\,Zm) be the space of paths which are right continuous and have left 
limits for some arbitrary fixed time T > 0. For a probability measure p'^ on Z^, we 
denote by the distribution on Z)([0, T],Zn) of the continuous time Markov chain 
rfit) generated by N'L^ with the initial measure f/'. 

The main theorem is stated as follows: 

Theorem 3.1. Let po : T' — > [0, 1] be a measurable function. If a sequence of proba- 
bility measures fi^ on Z^ satisfies that 



lim yu'^ 



-^yj(ON(x))r]^(0)- f J(u)po(u)du 

\Vn\ Jt^ 



> S 



= 0, 
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for every S > and for every continuous functions 7 : T'' — > R, then for every f > 0, 



lim I 

N—><x> 



y 7(<D^(x))77^'(f)- r J(u)p(t,u)du 

\Vn\ Jt^ 



>6 



= 0, 



for every 6 > Q and for every continuous functions 7 : T'' ^ R, where p(t,u) is the 
unique weak solution of the following quasi-linear parabolic equation: 

d 1 

-p = VDVp - 2rrn Z " p)'^v(.)//) , p(o, ■) = po(-). (3.1) 



Here we define Vv(p) := Yfi=\ Vi(e)(d/dxi) for e e Eq. 
We give examples corresponding to ones in Section IZ2l 
Examples 

0. The one dimensional standard lattice. 

For the embedding in Example Oa., we recover the equation in Theorem 3.1 in 

m-. 



d 



1 



d 



dt 



dH\ 



dx I 



For the embedding in Example Ob., we have the following equation: 



d 



5 



5 d 



dt^ Adx^^ 2dx 



p(l-p) 



dH_ 

dx 



1. The square lattice. 

For the square lattice and its embedding in Example la., we have the following 
equation: 



d 



1 / 



dt' 



2 \dx^ dy^ 



d 



dx 



dH\ d 



dx j dx 



dH 



dy 



For the square lattice and its embedding in Example lb., we have the following 
equation: 



A -11^ IJiL- IJL 

dt^ ~ \8 dx^ 4 dxdy "^4 dydx ^ 2 5/ 

5 d I dH\ I d ( dH\ I d ( dH\ d ( dH\ 

--4Yxr'-p^^j-2Mp^'-^^i;j-2Tyr^ 

2. The hexagonal lattice, the Kagome lattice. 

For the hexagonal lattice, the Kagome lattice and their embeddings in Example 
2. and 3., we have the following same equation: 

d 3/52 52 \ 3 5 / ^ dH\ 3 d dH\ 

dtf'^Ad^^^dfr'^d-^ r Ad-yr'-f"^^}- 
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4 Replacement theorem 



In this section, we fomulate the replacement theorem and give its proof. The replace- 
ment theorem is given by the form of super exponential estimate and follows from the 
one-block estimate and the two blocks estimate. 

4.1 Local function bundles 

For our purpose, we introduce local function bundles which describe the local interac- 
tions of particles and the two types of local averages for local function bundles. 

Definition 4.1. A local function bundle f on V xZ is a function f : V x Z — > R, which 
satisfies that for each z G V there exists r > such that f^:Z^R. depends only on 
[rix I d{x, z) < r]. Here d is the graph distance in X. We say that a local function bundle 
f : VxZ Ris T -periodic if it holds that fa-xic^) = fx{n)for any cr eF, x € V and 



Here we give examples of F-periodic local function bundles on V x Z. We use the 
first one and the third one later. 
Examples 

• If we define / : V x Z ^ R by for x e V and ?? e Z 

fxin) = Vx, 

then / is a F-periodic local function bundle on V x Z. 

• If we define / : V x Z ^ R by for x g V and 77 g Z 



then / is a F-periodic local function bundle on V x Z. 

Note that a F-periodic local function bundle / : V X Z -> R induces a map / : 
Vjv X Zat — > R for large enough in the natural way. To abuse the notation, we indicate 
the induced map by the same character /. 

First, for /? > 0, we define the R-ha\\ by 



7/ G Z. 




then / is a F-periodic local function bundle on V X Z. 



• Fix e G £0. If we define /'^^ : V xZ ^ Rby 




VocreVta-e there exists (unique) cr g F such that x = oae, 
otherwise. 



o-€Y,\cr\<R 
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Here | ■ | is the word metric appearing in Section |6T| We regard that B{Dxo,R) c 
via the covering map when is large enough for R. For a local function bundle / : 
y X Z ^ R, we define the local average of / on blocks {crDj;„)o-(=r by for x eV, 



where [x] is a unique element cr e F such that x € crDvo and denotes the 

cardinality of the set, which is equal to |B(Dj„,/?)|. Note that f^R- f^^n for every 
X e Djc„. As a special case, we define for 77 e Z and x eV, 

Second, we define the local average of / on each F-orbit, {crx}a-er by for x e V, 

Note that /. ^ and /,« are F-periodic when / is F-periodic. 

If a local function bundle / is F-periodic and is large enough, then /. ^, / « induce 
the functions on in the natural way. To abuse the notation, we indicate the induced 
maps by the same characters /. ^, f. R. 



4.2 Super exponential estimate 

For a local function bundle / : V x Z — > R, for jc € V, let us define {fx){p) '■- Ey^ [fx], 
the expectation with respect to the Bernoulli measure Vp. 

The following estimate allows us to replace the local averages of the local function 
bundle by the global averages of the empirical density. We call the following theorem 
the replacement theorem. We prove it in the form of the super exponential estimate. 

Theorem 4.1. (Super exponential estimate of the replacement theorem) 

Fix T > 0. For any Y -periodic local function bundles f : V X Z ^ R, for every 
x e Dx„ and for every 6 > 0, it holds that 

lim limsuplimsup-^logP^j-^ f Vx,N,e,K{ri(t))dt > 6 
e^o w-^oo |i n\ \ |i n\ Jo 

where 

Vx,N,,,k(tI) := X \fz.'c,K(v) - {fx)(Vaxo.,N)\ ■ 

Note that for every x e 0,^, t]^xo,r = '7o;.v,« > 0)- 

We denote by the corresponding distribution on D{[Q,T],Zff) of continuous 
time Markov chain rf^it) generated by N^Lf^ with the initial measure fi^. Furthermore, 
we denote by P^' the corresponding distribution on D([0, T],Zn) of continuous time 
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Markov chain rf^it) generated by N^L^ with the initial measure v^^j' i-^ ' ™ equilib- 
rium measure. We denote by the expectation with respect to P^, by En the one with 
respect to Pjv and by E)j^ the one with respect to P^. For a probability measure ju on 
some probability space, we also denote by the expectation with respect to ji. 

By the following proposition, we reduce the super exponential estimate for P^ to 
the one forPjy. 

Proposition 4.1. There exists a constant C(H, T) > such that 

-JL <sxpCiH,T)\TN\. 

Proof. To simplify the notation, put Hit, x) - H{t, <Siv(-^)) for x G V^. To calculate the 
Radon-Nikodym derivative: 



= exp 



exp 



xeVn xeVf/ 

Jo 

XT ( O TT 



exp 



By the inequaUty \e^-l-z- (l/2)z^| < (l/6)|2pel^l for z e R, we have that 

^ y^H(t,te)-H(t,oe) _ 

-AT^ ^ (Hit,te)-Hit,oe))-^N^ Y (Hitje) - Hit,oe)f\ 
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and thus, since (i> is harmonic, 

^2 ^ ^^Hit,te)-Hit.oe) _ jj 

1 d^H 1 

Furthermore, 

j^2^^H(,,,e)-H(,,oe) ^ ^H(,.oe)-H(Ue) _ 2) ^ (Vv(,)i/(f, Oe))' + O^. 

Hence, there exists a constant C(//, T) > depending only on H and T such that 
d¥^/d¥jii is bounded from above by exp C(H, T)\rf^\. It completes the proof. □ 

The super exponential estimate for induces the one for since it holds that 
for any Borel sets A c D([0, T],Zn), ¥"{A) < (expC(i/, r)|riv|)P/v(A). Furthermore, 
it is enough to prove the super exponential estimate for since for any Borel sets 
A c D([0, T],Zn), Pn(A) < l^^^^F'j^iA). 

Let 'P{Zn),'P{Z) be the spaces of probability measures on Zn,Z, respectively. De- 
fine := ®v„v[p, V :- ®vv\/2 the (l/2)-Bemoulli measure on Z^, Z, respectively. 
Here we introduce a functional on 'P{Zn), which is the Dirichlet form for the density 
function. 

Definition 4.2. For /i e !P(Zai), put the density (p — dfi/dv^. The Dirichlet form of 
is defined by 

In(m) = - [ ^^^LN^/^dv^. 

Note that hifi) = (1/4) /^^ ZeeE.i^. yl^fdv"^ > 0- 
Remark. 

The functional 1^ is also called the /-functional in the different literatures. 

Let us define the subset of the space of probability measures on Z^ by for C > 0, 

'Pn,c '■- |/" £ 'P{'^n) I /" is FAT-invariant and lN{pt) < C^^^j . 

The proof of the super exponential estimate for P^' is reduced to the following: 
Tlieorem 4.2. For every C > and every x e Dj^, 

lim HmsupHmsup sup \fx,K - {fx)(j]^a f^) = 0. 
First, we prove Theoren j4. ll bv using There m l4~2l 



Proof of Theorer mM Fix any T > 0. By the above argument, it is enough to show that 
for any F-periodic local function bundles / : V x Z ^ R, for every x e D^^ and for 
every 5 > 0, 

lim lim sup lim sup log P^M f V_,^N,e,KiMt))dt > 5 

'^^'^ e^Q N^oo \i-N\ \\In\Jq 
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where 

V^,N,,,Kiri) := ^ \l<TX,K{Tf} - {fx}('iiaxo,eN)\ ■ 

By the Chebychev inequaUty, for every a > and every 6 > 0, 



1 

1 n\ Jo 



T 





Now, the operator N^Lff + aVx,N,e,K acting on L^iZt^, v^) is self-adjoint for all a > 
and all x e Dj,,. Let AxM,£,K(fl) be the largest eigenvalue of N'^Lm + ciV_^^M,e,K- By using 
the Feynman-Kac formula, 



Jo 



exp a V,,^N,6,Kdt < exp 7/1^,^/,^,;^ (a). 



Therefore, it is suffice to show that for every a > 0, 



1 

lim limsupHmsup -r—^x.N,£.K{a) - 0, (4.1) 

K^'" e^O N^ca \^n\ 



since ( I4.1I ) implies that 



1 



1 

\^n\ Jo 



lim lim sup lim sup -pp-;- log j j V x^i^^^Kiri{t))dt > 6\< -ad, 



and we obtain the theorem by taking a to the infinity. 

By the variational principle, the largest eigenvalue ■Ijc^n.€,k(<^) is represented by the 
following: 

yix,N,<!,K(a) = sup \a \ Vx,N,6,Kdfi - N^In(m) 

See fT\ Appendix 3 for more details. 

Denote the average of ^ by F^-action by Ji :- (l/IFivl) Xo-gFa. ° £!■ Then Jl is 
FAf-invariant so that 



1 r 1 |~ 

TTT- V,,N,e,Kdn = — [Vx,N.>;,k] = % \f,,K ' {fx 

|i n\ Jzn m ' 



The functional In(-) is also FAr-invariant, i.e., Iff(jj o o;) - I/^dS) for jj e ViZt^), 
a € Fiv. Thus, it is suffice to consider FAf-invariant measures /j. to estimate the largest 
eigenvalue /ij:,iv,f,A:(fl). Furthermore, it is suffice to consider the case where yu satisfies 
Vx,N,c.KdiJ > N^In(ij)- There exists a constant C(f) > depending on / such that 
Vx,N.e,K < C(/)|Fai|, thus we reduce fi e ViZ^) to every F/vi-invariant measure satisfy- 
ing that for every C > 0, /wCju) < C\Vn\IN^- This shows that we can reduce to 'Pn,c for 
every C > 0, and thus it is enough to show for every C > and every x e D j„, 



lim lim sup lim sup sup U^^k - ifxXrixo ew) = 0, 

f^O //ePwr ' 

to obtain ( 14. It . This follows from Theorem [4.2| It completes the proof. 
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4.3 The one-block estimate 



In this section, we prove the one -block estimate. We regard a probabiHty measure yU 
on Z^r as one on Z by periodic extension. Let pr^r : V Vjv be the covering map by 
AT- action, and define the periodic inclusion ip^^ : '—^ Zby iiper^)z '■= Tlprn(z),'n £ 
Zm,z e V. We identify // on Z^ with its push forward by i*^. On the other hand, we 
identify an A^F-invariant probabihty measure on Z with a probabiltiy measure on Zjy. 

First, for a finite subgraph A = (Va, Ea) of X, we define the restricted state space 
Za := {0, 1}*''^ and the (1 /2)-Bemoulli measure by :— (giy^vj^j ^a- Let us define 
the operator acting on l}{Z,v) by := (1/2) 2c££a ^c- For n e P{Z), ;u|a stands 
for the restriction of fi on Za and 0a := d/^lA/dv^ its density. We also define the 
corresponding Dirichlet form of V0A by 

JZa 

For large enough A', we regard A as a subgraph of X^r by taking a suitable fundamental 
domain in V for AT-action. 

By the convexity of the Dirichlet form. 




by putting := d^/dv'^. 

The one-block estimate is stated as follows: 

Theorem 4.3. (The one-block estimate.) For every Y-periodic local function bundles, 
f : y X Z — > R, every x e D^^ and every C > 0, 

lim lim sup sup E^^ \%,k - {fx) (^;co,/i:) = 0- 

Proof. For any probability measure € f Ar,c, we apply the above argument by setting 
A as Va := B{Dx„,K) and Et^ := [e & E \ o{e),t{e) e Va}- Since ju^v and are F^y- 
invariant, 

< \{a e I \a\ < K]\ ■ -^In(mn). 

\^n\ 

Since e P^^c and |Vivl/|FA,| = IVqI, it holds that 

We note that fiZ) is compact with respect to the weak topology, and thus {//jv) c 
PiZ) has a subsequence which convergences to some fi in PiZ). Let ^ be the set of 
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all limit points of {/un} in 'P(Z). By the above argument, /^(yu) - for all /j. e 
Therefore, since /°(//) = (1/4) ^^g^^ I„,ezJ'^e yjl^kiri))^ = for every e ^ by the 
definition of we obtain that //|a(?7^) = /u|a(?7) for every e e and every 77 e Za. 
This shows that random variables 771, x e V are exchangeable under By the de Finetti 
theorem, there exists a probability measure A on [0, 1] such that yu = VpA{dp). Since 

Um sup sup \Z,K - (fx) {rix,,K)\ ^ sup E^, \f,,K - {fx) (^x„,if)| . 

it is enough to show that 

Umsup sup Ey^ \f,,K - (fx) ('?t„,A:) = 0, 

K^oo p£[0,l] ' ^ /| 

for every F-invariant local function bundles /. 

By the definition of the F-periodic local function bundle, there exists a constant 
L > such that for every x e V, fx '■ Z ^ R depends on at most {rj, \ z e L)}. 
Therefore we obtain that there exists a constant C(f) > depending only on / such 
that ^ 

Ev, [fx,K - Ev, [fx,K]f < Cif) ■ |g ^ as ^ cx,. 

Note that {fx)(p) - Ey^^ [/x.at] for every x € D^^ since the Bernoulli measure Vp 
is FAT-invariant. In addition, we also obtain that there exists a constant C > not 
depending on p, 

|2 C 

Ey, |?7xo,jf - p| ^ ^ as ^ 00. 

Finally, since {fx){p) is a polynomial with respect to p, in particular, uniformly contin- 
uous on [0, 1], we obtain that 

lim sup Ev,, \f,,K - </t>(^i-„ k)\ = for every x e £)^-„. 

This concludes the theorem 

lim lim sup sup E^„ \f,^K - {fx)(r]xo k)\^^ for every x e D^^. 

□ 

4.4 The two-blocks estimate 

In this section, we prove the two-blocks estimate. We identify a probability measure 
on Za with its periodic extension on Z in the same manner as Sectior l43] 

Theorem 4.4. (The two-blocks estimate) For every C > 0, 

lim lim sup lim sup lim sup sup sup E^fjljli^ k "'Ho-xo k\ - ^■ 

^^'^ e^O A-*oo a-ers.t.L<\a-\<t;N ^^HePf,.c 
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Let us denote by !P(Z x Z) the space of probability measures on Z x Z. We define 
the map for creF, (t:Z—>ZxZ by &(ri) := (//, err]), rj e Z. For fi € !P(Z), we define 
the push forward of yU by o" via &ii :- jio o-"' e !P(Z x Z). 

Let us denote by c !P(Z x Z) the set of all limit points of [d-jiN \ L < \a\ < 
eN,fiN e Pn.c) in PiZ xZ) as N ^ oo and by ^° c P(Z x Z) the set of all limit points 
of y[° as L — > oo. We put (xq, x'q) e V x V, where x'q stands for a copy of xq. Then, it 
holds that 

I |'?x„,if - ^^jco.kI t^idr]) = I k^^ i. - 1]'^, J (d-fi){dr]dT]'). 

Jz JzxZ ' ' 

Note that 

UmsupUmsup sup sup I \Tj^^i^-r][,J(&f/)(dr]dT]') 

L-^oo A'-<(XD o-ers.t.L<\a-\<eN fiePn.c JzxZ ' ' ' ' °' ' 

^ sup I |77,^_^ - 7;;, ^L(t/77t/77'). 
Jzxz ' ' 

We introduce two types of generators acting on L^{ZxZ, v®v) and the correspond- 
ing Dirichlet forms. The first one is used for treating two different states at the same 
time independently. The second one is used for treating exchanges of particles between 
two different states. 

As in Sectior |431 we define a subgraph A = (Va,Ea) of X by setting Va '■- 
B(Djc„, K), Ea := {e e E \ o(e), t(e) e Va) and the operator acting on L^(Z, v) by 

For yu e !P(Z x Z), we denote by ;u|axa the restriction of yu on Za x Za. Define the 
Dirichlet form of V^axa by 

JZt^xZf, 

Let us introduce the notation which describes exchanges of states for {r], rj') e ZxZ. 
For ix,y) e VxV, (rj, rj'Y^'^^ e ZxZ is the configuration obtained by exchanging values 
77, and T][., i.e., (77, t/')*""'-*' '■- ((zy('z'\z,z')evxv is defined by setting 

yr], z + X, 

Moreover, for F e l}(Z xZ,v ® v), we define tt^jF ((tj,1]')) :- F ((//, t/')*"''^'') - 
F ((t], rj')). We define the operator acting on l}{Z x Z, v (gi v) by 

I- o 
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The corresponding Dirichlet form of V0axa is defined by 

We prove two lemmas needed later. The first one is easy to show, so we omit the 
proof. 

Lemma 4.1. There exist constants ci,c\,C2 > such that for all o" G F, 

c[\cr\ - C2< d{xQ, axo) < ci\a\, 
where d is the graph distance ofX. 

For X, y G and for t] g Z^t, 77^*'^' is the configuration obtained by exchanging two 
values T]jc and i.e., 

% z-x 
Tlx z = y 
T]^ Otherwise, 



and moreover for x,y € Vn, we define the operator nx,yF(j]) := F^t/*^-^') - ^(77) for 
F e L^(Zn, v^). These notations also indicates ones for Z. 
The second one is the following: 

Lemma 4.2. For every Tt^ -periodic functions F g I^iZf^, v^) and every a g F^t, 

r {7:^,a.,Ffdy^ <Ad{xo,axQfYj f i^^Pfdy^ . 

Proof. Forxo,crxo G W, there exists a path c = (ei, ... ,eO such that xq = o(e\),t(ei) = 
oyco and d{x(), o^xq) = Define a sequence of edges c" = (/i ,...,/;, , ... , /iz+i) by 
setting c" = (ei, ■ . ■ , e;, . . . , ei). For 77 G Za?, let us define //(O) := 77, 770) := ?7;^'_i), 1 ^ 
J < 2/ - 1, inductively. We note that rf^-°''^ = t](21-i)- Then we have that 

'21-1 21-1 
fco.cr.o^'C'?))' = (F irfo''^'^) - F iri)f = ^ imi-i)) < (2/ - 1) ^ {^A ('70-i)))' • 

V i=l ) !=1 

If f is Fjv-periodic, then for each i = 1, . . . , 2/ - 1, 

r {nf,F(rJ^i_,))fd^^<Y,\ i^eFfdv^. 

Therefore 

•JZk JZf/ JZn 

It completes the proof. □ 
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Let us define the subset of ViZ x Z) by for a constant C > 0, 

^e,c := {/^ e P(Z X Z) I I^^^itx) = 0, < Ce^} . 

Then we have the following lemma. 
Lemma 4.3. There exists a constant C > such that 

c ^. 

Proof. As in Sectio H431 we define a subgraph A := (Va,£'a) of X by setting Va :- 
B{D^„,K),Ef^ {e e £■ I o(e),t(e) e Va) for K. We take large enough L,N for the 
diameter of A, K, so that VAno-VA = for L < |cr| and VaUo-Va c Xa, for L < |cr| < eN 
by taking a suitable fundamental domain in V by AT-action. Take jj e 'Pn,c, ct e F. 
For /?') e Za X Za, if we define ^ e ZauctA by T/Ia = and ct^^tjIa - t]' then 
(a-yy)|AxA('7, ??') = A'Iauo-aC^- Let us consider the operator acting on L^(Z, v) by 



1 

-^auo-a - 2 



AUo-A 



For pi e !Piv,c, we denote the density of /ilAutj-A by (pi\}ja-K '■- dfi\AU(TA/dv . Then 

^Axa('^/^) = '^AUcrA(/^) = - I V^AU^^AUtrA yJ'pAUcrAdv^'^"'^ . 

For any // e !PAf,c5 we put := dp/dv^ . By the convexity of the Dirichlet form and the 
FAT-invariance of (p, 



Since Im(m) < C\Vn\/N^ 



JO \Ea\Wn\ C \Ea\,,.,C 

W-/^)^— ^p-— as A. ^ CO. 

Therefore = for any /Ji e Furthermore /^xaO^) = for any ju e yi° 

by the continuity of the functional ^axa- 

For fi e !Pai,c> by the convexity of the Dirichlet form. 

By the FAr-invariance of // e ^Pat.c and by Lemm d4~2l 
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By Lemm d4~T] and the definition of yU € Vn,c^ for all cr e F such that L <\a-\ < eN, 

(.x^) < 4c>i^ ■ ^ ■ < 

By setting C = 16cjC|yo| and the continuity of the for every Ji e we have 
that '® ^ f^C. It concludes that y[° c ° 

Let us prove Theorer d4~4l 
Proof of Theoren W^ Denote by JIq the set of all limit points of J^^q as e — > 0. For 
every Jlo e yio, it holds that /^xa^o^ ~ ^ ^axa°'^o^ = by the continuity of the 
functionals /^^^ and /J^'"^''^ These show that for any x,y e Va x Va, T^x.y (J^oIaxa) = 
and thus for any x,y eV^ and for any 77, 77' € Z, ^uoIaxa {{rf'''^\ rj')) = i«olAxA (('?, ??')), 

MoIaxa ((??, 77'^'''^^)) = A^oIaxa ((??, ??')) and JuoIaxa ((??, t?')*'*'''^) = I^oIaxa ((??, ??')), i-e., /To is 
exchangeable on Z x Z. By the de Finetti theorem there exists a probability measure A 
on [0, 1] such that/7o = jj Vp ® VpA(dp). 

^ I i2 

As in the proof of Theorem [4.3l lim/f^oo sup^^jg i] I'? a: ~ p| - 0' therefore, by 
the triangular inequality, 

sup p7,.^_^ - ^ < sup Ey^^e>v„ p7x„,if - Vx'kI 

-poema ' ' pe[0,l] ' ' 

< 2 sup Ey^ It;^^ a: ~ p| ~* as /iT 00. 

pe[0,l] 

Finally, by LemmjOl ^° c ^ for some C > and thus, 

lim sup lim sup lim sup sup sup ^ - rj^^-, J 

e->0 L->oo N-icx, a-ers.t.L<\a-\<t;N fiep^.c ° 

< Hm sup sup E^ - Tj'^, J < sup E^„ 1^^^^ - Ij'^, J -> as ^ 00. 

This completes the theorem 

lim lim sup lim sup lim sup sup sup E^ mj^^ - rj^^j A-Q. 

e^O L^oo W^tx) iT£Ys.l.L<\tT\<eN iiern.c ° 

□ 

4.5 The proof of TheoreitM 

Let us prove Theorerr f4.2| bv using the one-block estimate Theorerr |4.3| and the two- 
blocks estimate Theorerr f4.4l 
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Proof of Theoren {4~2\ First, we note that there exist positive constants Cq,C\ and C2 
such that for any rj eZ, 



1 



\^L<\cr\<eNO-D 



Z 



id K'' [(eN)'' - {eN - 1)"') K'' ((L + 1^ - U') 

^ Co 7 + C 1 ; h Co ; , 

and thus there exists a constant C(e, L,K) > depending on e, L and TiT such that for 
any rj e Z, 



1 



'7.10, eA' 



C(e,L, 



uniformly. Then, since ju is F-invariant as a probability measure on Z, 



1 



1 



^L<\a-\<6NCrDxA 



Z 



UL<|o-|<eWO"D 



Z 



C(6, L,ji:) 



^1- - I C(e,L,K) 
< sup sup \t],^ i. - 77^,^^| + , 

a-£r,L<\a-\<eN fiepN.c ■'^ 

where the last inequality comes from the fact that for any z e o-D^^ it holds that rj^.^ = 

\.K- 

Applying Theorerr f4.4l we have that 

lim lim sup lim sup sup E^ \Tj^^ ^ - rj^^ - 0. 

For every F-periodic local function bundles /, ) is uniformly continuous on [0, 1]. 
Therefore, 

lim lim sup Um sup sup E^ </v)fe„,;r) - </t) (^x„,<,Ar) = 0. 

Furthermore, applying the one-block estimate Theoren j4.3J for every x e D^^, 
lim lim sup lim sup sup E^ - </J (^,^ = 0. 



It completes the proof of Theorer rf4.2l 
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5 The proof of TheoremlXTI 

In this section, we prove Theoren f3.1l 

Let : Xn — » T'' be the A^-scaling map. We define the empirical density by 

^N(t,dfi) tJ^ y 77^'(0<5o„w(af/i), 

where 6^ is the delta measure at z e T''. The empirical density is the measure val- 
ued process. We denote by C''^([0, T] x T'^) the space of continuous functions with 
continuous derivatives in [0, T] and twice continuous derivatives in T'. For every 
/(■) e C''2([0, T] X T''), we define 



To abuse the notation, we denote the inner product in L^{T^, v^) by 

{F, G) = r FGdv^ for F, G e L^(T', lA). 
Jt'' 

Let us define the process as follows: 

Jo 

where ^^(f) {id/dt)J„^Nit)) + N^L^{J„^N(t)) and 

Jo 

where 

Here Mff{t),NN(t) are martingales with respect to the filtration {;F,),>o, where := 
cr{77'^(s) I < < f! and it holds that 



Jo 



Then we have the following lemma by applying the Doob inequality. 
Lemma 5.1. 

^ 0. 



hm E?, 



sup |MAr(f)|^ 
0<r<r 
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Proof. For J e C''^([0, T] x T'') and for each e € E^v, we have that 



X€Vn 



IV, 

By the regularity of J and the compactness of T'', we note that there exists a con- 
stant C( J) depending only on J, such that uniformly, 



\Js(^N(te))-Js(^N(oe))\< 



N ' 



Thus 



AN{s)< — YjC"{.e,r],s) 



e€En 



1 C(J) 



\Vn\ N 



, , C(JY C" 



where C is a constant such that c"{-, ■, ■) < C and C" = C'C(7)2|£:ol/|Vo|. Then we 
obtain that 



Jo 



C"T 



as A' ^ oo. 



JO \Vn\ 
Applying the Doob inequahty for the right continuous martingale {M;v(0}«> 



sup \MNiT)\' 

0<t<T 



<4E^\M!,iT)\\ 



we conclude that lim^v^oo E^ l^supo^f^j- |Miv(/)P] =0. □ 
5.1 Relative compactness of a sequence of probability measures 

We denote by := A1(T"') the space of nonnegative Borel measures with the total 
measure less than or equal to one on T'', endowed with the weak topology. Since T'' 
is a compact metric space, the space of continuous functions C{T^) with the supre- 
mum norm is separable. Fix a dense countable subset {Jk]'^-Q of C{T^), then the weak 
topology of M is given by the distance dMi-, •) by 



dM(M,M )'■= y,^- T— 177 — — 77 — 7 

^2* \ + \{Jk,n)-{Jk,ti 



for e At where {Jk,iJ) '■= Jkdfi. We note that At with the weak topology is 
compact. 

Define the space of paths in At by 



D([0, J], At) := {^. : [0, J] ^ At I ^ is right continuous with left limits.} , 
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equipped with the Skorohod topology. For a given process '■ [0, T] ^ M such that 
i^N e D([0, T],M)) = 1, we denote by the distribution of on D([0, T],M). 
Then we show that the sequence {Q^}n has a subsequential Hmit. The following 

proposition gives a sufficient condition for this. See Q Section 4, Theorem 1.3 for the 

proof. 

Proposition 5.1. If for every J^, k - 0,1, . . ., and every 6 > 0, 

= 0, 



lim lim sup 

7~>0 AT^co 



sup \{Jk,Mt)) - {Jk,Ms))\ > S 

Wi-s\<y 



then there exists a subsequence {Qj^^J^Q and a probability measure on D{[Q, T~\,A\) 
such that weakly converges to as A: — > oo. 

The next proposition claims that each subsequential limit in Propositior lS.ll is 
absolutely continuous with respect to Lebesgue measure on the torus for each time t, 
and its density has the value in [0, 1] a.e. The proof is the same as in |7| Section 4, 
pp.57, so we omit the proof. 

Proposition 5.2. All limit points o/{(2j^)iv ore concentrated on trajectories of ab- 
solutely continuous measures with respect to the Lebesgue measure for each time t, i.e., 
there exists a Borel set W c D([0, T],M) such that Q^iW) - 1 and for every ^. e W 
and every t € [0, T], is absolutely continuous with respect to the Lebesgue measure 
du. Moreover, the density p(t, u) :— d^,/du satisfies that < p(t, u) < 1, du-a.e. 

To simplify the notation, we put J{t,x) :- J{t,<^j^{x)),H(t,x) := //(r, <E)a-(ji:)) for 
J, He C^'HlO, T],T''), respectively. 
Then we have: 



9 H ^ 

N-'L^iJ,, ^N(t)) = > ^ exp [Hit, te) - H(t, oe)] ■ i-T]oeT],e + T]oe)(J(t, te) - Jit, oe)) 

= :;T^ X y! ii^^P Wt, te) - H(t, oe)) - 1) ■ (J(t, te) - J(t, oe)) rj^e + (J(t, te) - J{t, oe)) rjoe] 



^ {exp (H{t, te) - H{t, oe)) - exp (H(t, oe) - H{t, te))} ■ T]oer]te{J(t, te) - J(t, 



oe)). 



For e e £, we denote the directional derivative along v(e) by 

i=l ' 

Applying the inequality |e~ - 1 - z| < (l/2)|zpe'^' forz € R, by the regularity of H 
and by the compactness T'', there exists a constant C > not depending on each point 
of [0, T] X T'' such that for every and for every e e E^, 

\N\exp {Hit, te) - Hit, oe)) - 1) - Vv(,)//(f, oe)\ < |, 
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|A^{exp (H(t, oe) - H(t, te)) - 1) + Vv(,)//(f, oe)\ < 



C 



\N{J(t,te) - J{t,oe)} - Vv(e)7(f, oe)| < — . 
By the convergence of the combinatorial Laplacian in Sectior lZSl we have that 



^ (J(t,te)- J(t,oe))^ - ^ Y^-—-(t,x)vi(e)vj{e) + OM 



(?e£A',A- 



Nj l.J=l ■' 



Hence, 



N^L^{J„^N(t)) = t-Jt- y Vv(,)//(f,oe) ■ Vv(,)y(f,oe)77„ 



eeE^ 



-(t,x)vi(e)vj(e)j] 

oe 



1 

dXidx: 

' ' xeV^ eeEfij- i,j=l ■' 

~ ^^7777 y! '^x(e)H{t, oe) ■ Vv(e)7(f, oe)rioeT],e +on. (5.1) 



We replace Pj^ by regarding the empirical density as the measure on Z3([0, T] , A1). 
Let us prove the following lemma: 

Lemma 5.2. For every J e C(T'') and for every 6 > 0, 



lim lim sup 2^ 



sup <7,^^(f)>-a,^ivW> >5 

|/-.s|<y 



0. 



Proof. For every continuous functions 7 : T'' — > R, it holds that 

Since by (15. It there exists a constant C such that for large enough A^, bN(t) - {(d/dt)J,, ^^(1))+ 
N^L^{Ji, ^Nit)} < C uniformly, we obtain that by the Chebychev inequality and by the 
triangular inequality, for every 6 > Q, for every y > and for large enough A^, 



Then by Lemm dSTl 



sup \{J,^N{t)}-{J,^N{s)}\ >6 

\t-s\<y 



Cy + 2 sup |MA-(f)| 
o<t<r 



lim sup Q" 



Therefore for every 6 > 0, 



sup \{J,^N(t)}-{J,Ms)}\ >S 

\t-s\<y 



<cl. 



lim lim sup 



sup \{J,^N(t)}-{J,^N(s)}\ >6 

l'--5|<y 



= 0. 



It completes the proof. 
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5.2 Application of replacement theorem 

We have the following estimate: There exists a constant C{H, J) > depending only 
on H and J such that for every e e and for every a^, a' e r^, 

<C(i/,7)|k'-^||ii. 



Thus, by the regularity of H and J, putting GNipe) :- 'Vx(e)H(t, oe) ■ 'Vx(e)J(t, oe) for 
e e En, 

V' V' \Gisi{ooy) - GNoaeA ^ C— for some constant C > 0. 

Here we regard G(-) as a local function bundle independent of states and denote by 
Gffoo-e.K the local average. By the uniform continuity of the twice derivative of J, 
putting 

F^ix) 2^ 2^ g^ g^ (^-^^)v,(e)vXe) for;c e Vn, 



it holds that ^ 



Here we also regard Fm( ) as a local function bundle. 
By the above argument, we obtain that 



Here Z**^* is the local function bundle appearing in the third example in Sectioij4T| 
and /('"^ ^ its local average. 

By applying Theorerr f4.1l and by the continuity of (d^ / dxjdx i)J , Vv(e)H and Vv(e)/ 
on the compact space [0, T] x T"', it holds that for every t e [0, T] and for every <5 > 0, 

lim limsuplimsup-^logP^f f Itt^ V, X ^'^(Z-^) fe.v./f - ^^^.^.^w} xj) = - 

e^O N^^ ir^vl \ Jo I \Vn\ ^ - - ' I / 



lim Um suplim sup log f I - j- V V GA.(oo:e) {770^,,^ - ?7„„,^^} \ds 



> 5 = -00, 
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and 

lim limsuplimsup -L logP^( f iJ- V V GNiogy) [f'\^e^K ' (^o^^.^iv)^} k > ^ 

Here we use {f!f^}(p) - for every x e in the third estimate above. By the 
triangular inequality, it holds that for every t G [0, T] and for every ^ > 0, 

lim sup lim sup f \n'-L^{J„Ms)) 

e^O N^oo \ Jo ' 



2\V. 



1 d^J 



TiTm Z Z ^mH(s,ooy;) ■ V^(e)J{s,oae) ■ (n^^^^^^^fUs > = 0. 
Applying the convergence of the combinatorial Laplacian in Section [231 we have 



that 



IT^ Z Z Z y ^^(^o:-^>;(e)v>(e)= y VDV7(f,o:jco) + Oiv. 
Wn\ ^ dxidxj \Yn\ 



Recall that 



<7,,^M(f)) - </o,^^(0)) = j| {<5,7,,f^v(*)) + N^L'^{Js,^n{s))] ds + MM. 
By Lemm jSTTj and by the Chebychev inequality, for every 6 > Q, 

e« I sup \MN(t)\ >6]< (1/(J)E^ ( sup |M^(f)ll ^0 asN^oc. 
V)<t<T I \0<r<r / 

Furthermore, by the triangular inequality, we have that for every 6 > and for every 
t e [0, T], 



lim sup lim sup 



<7„^m(0)-</o,^w(0))- I \{dsJs,^N{s)) 



" TlFT Z Z ^y(e)H{t,ooy) ■ Vv(.)7(f,oo:e) ■ (^^^„,,a,)^ > 5) = 0. 
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By Lemm j6^ we replace rj^j.^^ by {^N,X<bN(s:.->:o),<;) ^"'^ 'he summation for cr e F^v 
by the integral. Since by Lemm d572l and by Propositior j5.ll the sequence {2^1 is rel- 
atively compact in the weak topology, for a limit point there exists a subsequence 
{QnJ weakly converging to Q^. By Propositior |5.2| the empirical density concen- 
trates on an absolutely continuous trajectory p,djj as — > oo. By the assumption of 
Theoren im we replace {Jo,^NiO)) by {Jo,Po), and then we have that for every 6 > 
and for every t e [0, T], 



Um sup 2^ 



{Ji,Pt} - </o,Po) 



+ ^y. r '^Me)H(t,Z)-'7,(e)J(t,Z)-{p„Xz,Mdz)+ f VBV J(t, z) ■ {p„Xz,Mdz) 

A vol 5r -JT'i JT' 



>(5 =0. 



Vy(,)H{t,z) ■ Vy(,)J{t,z) ■ {p„Xz,,) Midz)>ds 



Here we replace v(e) for e e E'^ hy v(e) for e e £0 by the F-invariance of v(-)- 
By the Lebesgue dominated convergence theorem as e — > and by the triangular 
inequality, we have that for every 6 > and for every f e [0, T], 



{Jt,Pt) - <-/(), Po) 



r \{dsJs,Ps) + y <Vv(.)7 ■ Vvm//,a) + <VDV7,p,) 



This shows ^N^it) concentrates on p, which is a weak solution of the quasi-linear 
parabolic equation (13.1b . Furthermore, p has finite energy by Lemm dTTI By the 
uniqueness result of the weak solution Lemm dX2l in SectioijT] we conclude that the 
limit point of {Q^}n is unique and concentrates on pdi-i as goes to the infinity. 
That is, for every 6 > 0, 

lim Q^{dst{^N,pdfi)>6)^0, 

where dsk. is the Skorohod distance in D([0, r],A1). In particular, since T > is 
arbitrary, it follows that for every f > 0, for every (5 > and for every continuous 
functions J e C(T''), 



lim 



y y((&^W)77^'(f)- r mp{uu)ii(du) 



> 6 



= 0. 



It completes the proof of Theorem lTTI 
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6 Appendix;A 

6.1 Approximation by combinatorial metrics 

Take a Z-basis cri, cr^/ of F and identify F with Z''. Define the standard generator 
system of F by 5 - [cri ,.. ., crd, -(T\ , . . . , -cr^). We introduce the length function 
associated with 5 , | ■ | : F — > N by 



for cr e F. Then the map {(t,(t') € F x F i-> \(t - cr'\ e N induces the met- 
ric in F, which is called the word metric associated with S . Let crj , . . . , cr,/ e F^y 
be the image of a-\,...,a-4 by the natural homomorphism F — > Fat. Then S_ :- 
{cr\, . . ., (Td, -cTi, . . ., -cTd] generates Tn- The length function associated with 5_, | ■ | : 
Fjv ^ N is also defined in the same way. To abuse the notation, we denote the word 
metric in Tn associated with S_ by the same symbol | • |. 
We define an Z' -norm in F ® R = R'' by 



for X = {x\,...,Xd) e R'' and the distance d\ in R'' by (ii(x,x') :- ||x - x'||i for 
X, x' G R''. Denote by d\ the induced metric in T'' from di . 

Fix xq e V and a fundamental domain D^^ c V such that xq g D^^ and D^^ is 
connected in the following sense: For any x,y e Di„ there exist a path ei, . . . ,e/ in £ 
such that oe I - x,tei = y and oei, fei, . . . , oe/, fe/ are all in Dvo- This kind of set 
always exists if we take a spanning tree in Xq and its lift in X. To abuse the notation, we 
denote by xq e W, c Vn the images of xq, D.i-„ by the covering map, respectively. 
We also fix a fundamental domain c E which is identified with Eq. To abuse the 
notation, we denote by c E^ the image of E^ by the covering map. 

We define the map [■] : V — » F as follows: For x eV, there exists a unique element 
cr G F such that x G crD^^ since F acts on X freely. Define {x\ :- cr. Since there exists 
a constant Cq :- maXv^y ||<l>([x]xo) - <I>(x)||i by the F-perodicity, we have that 

\m[x]xQ) - <l)([z]xo)||i - 2Cq < mx) - <b{z)\h < \m[x]xQ) - <D([z]xo)||i + 2Co 

for any x,z G V. Furthermore, since ||0([x]jco) - 3)([z]xo)||i = \[x] - [z]\, we have that 

I W - [z]\ - 2Cq < mx) - 0(z)||i < \[x\ - [z]\ + 2Cq. 

As in Sectior jm suppose that we have an injective homomorphism i/r : F — > R'' 
such that 

il/(Y) - {/ kiUi I k; integers and m, g R'', mi, . . . , are linearly independents . 




d 



llxlli :=2n 



!=1 
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Take a fundamental parallelotope P :- |Ef=i <^iUi | < a,- < l| c F i8> R. In the similar 
way to the above, we define the map [•] : F igi R — » F as follows: For x e F gi R, there 
exists a unique element cr e F such that x e crP. Define [x] := <t. Since there exists a 
constant C\ := maxx^roR l|3*([x]xo) - x||i by the F-periodicity, we have that 

\m[x]xo) - a)([z]xo)||i - 2Ci < ||x - zlli < \m[x]xo) - a)([z]xo)||i + 2Ci 

for any x, z e F ® R. Furthermore, we have that 

ILxJ - LzJI - 2Ci < ||x - zlli < |[x] - [z]| + 2Ci. 

Let us define an e-ball in T'' of the center z e T'' by := {x e T"^ | di(x, z) < e}. 
We show the following: 

Lemma 6.1. There exists a constant C(e) depending only on e such that for any z G T'' 
and for any N >\, 



vol[Bl{e)) 



VOl{T'^) 



C{e) 
N ' 



Here vol{T) stands for the volume of a Borel set T and \U\ the cardinality of a set U. 

Proof For any z e T'',Jake a lift z e R'' then ||<l)([z]xo) - z\\i < Ci. For sufficiently 
small e > 0, take a lift B\(e) c R'' of B\{e) c T"'. Again, from the above argument, it 
holds that U|o-|<,/v-2CiCr[z]P c BliNe) c U|o-|<ea,+2c,o-[z]P. 

Note that vol (e)) = vol and thus |a^''vo/ - vol (U|<^|<eiva-[z]P)| < 

vol{P) {l'^ld\) [ieN + ICiY - (eN - 2CiY). 

Since |{o- e F | \cr\ < eN}\ = vol (U|^|<eA,o-[z]P) /vol (P) = |U|^|<eA,o-[z]Z);,J /\Vo\ and 
it concludes that there exists a constant C{e) > depending only on 

V0l[Bii6)) |U|^|<,;vO-[z]D;co| 



e such that 



vo/(T') 



Cje) 
N ' 



The cardinaUty of the set U\a-\<eNcDjCf^ is invariant under translation. It completes the 
proof. □ 

Let us define a measure // on by jU := (l/vo/('F'))<ix. hetXz,€ '■ T' — > R be a 
characteristic function defined by 



onT''. 



For the empirical density '■= (1/l^wl) 'LxeVn Vx^^nU) on T', 77 e Z^, then we have 
the following lemma. 
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Lemma 6.2. There exists a constant C{e) > depending only on e > 0, such that for 
any t] e and any z € Vn, 



where 



i^N,X^N{z),e) - 7- — — r Yj 



C(e) 

N ' 



Proof. Take a lift 2 e V of z e V^, and a lift £',^^^^(,,(6) c T ® R of ^ ^'^ ■ I" 

the similar way to the proof of Lemm dSTl we obtain that 



y cr[z]Z),„ c |x € y 



|o-!<eA'-2Co 



-cl)(x) <D(2) 



|<r|<eA'+2Co 



Furthermore, we take a lift 77 € Z of 77 e Z^^, then it holds that 2 o (aogb" (e) '7a 



A-ey,||iO(x)-i«(2)||,<E 



- y 



A-eUjA._2C|,<|ir|<6«+2Co o"[2]£*,-, 



The last term is bounded by (|y()|/|yAr|)(2''/£/!) ((eA^ + IC^if - {eN - 2Co)''), and thus 
there exists a constant Ci(e) depending only on e such that 



IVivl X \Vn\ X 



Ci(e) 
N 



ByLemmiOand|U|^l<^A,cr[(l/A?)(D(z)]£)^„| = |U|^|<,ivcrD,„| = |U|,^|<,ivo-[z]D.,„|, 



C(6) 



Wn 



where C2(e) is a constant depending only on e. Finally, 
1 



1 £K£)^C^ . 
< ^ 



where C3(e) is a constant depending only on e. It completes the proof. 
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7 Appendix;B 

7.1 Energy estimate 

In this section, we prove the following lemma. 

Lemma 7.1 (Energy estimate). Suppose that {Qn}n>\ is a sequence of probability 
measures on D{[Q,T],M). For any limit point Q* of {Qn]n>i, Q*-a.s. there exists 
a measurable function p(t, u) such that = Ptdn, p has 



for every J e C^'HLO, T] x T"*) and \<i<d. 

Proof. For fixed xo e Vn, we define a lattice Xf^ - (V;^, Z?^) whose vertex set V^i is the 
subset of Vat in the following: V'^ is the orbit of xq by Fjv, i.e., Fjvxo c Vn- Define £^ 
the set of oriented edges / such that / = (xo,o^xo) for some € {a;., -o\ \ I < i < d). 
Then F/v acts on naturally. A configuration on X/v induces the one on by 
restriction. We use the same symbol 77 for this restriction. For J G C''^([0, T\ x T'') and 
for J e {1, . . . , c?}, we define oF> : C\[0, T], R), 



—peL\[0,T]xT'), 



for I < i < d, and satisfies 





J{o f) if there exists g F^r such that / = a;(xo, cr,xo) 
otherwise, 



and 
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For any fi e P(Zj^), we have that 

/ I \^ 

^Zn JZn ^.ggs y y uv J 

{ I \^ 



< 



We use the Cauchy-Schwarz inequality in the last inequality. By using the argument in 
the proof of Lemm d4!2l there exists a constant C such that 



Since /a,(;u) = Y.eeE^ (^f y[dijjdv^^ dv^ and ZfeE'^i(^f(r]of-r],f)f < 2 Zaer^ Jigycaf, 
get 



7(^Xo)2 ■ ^j\6ClN(il). 



Note that Yjo-er^ Jigjio)^ < 2\Yn\ ■ WA^^i^jd-^ for large enough A^. Consider for a > the 
self-adjoint operator 

and suppose A^^(a) to be the largest eigenvalue of this operator By the variational 
formula 

A'-^(a) = sup 
ijenz. 



3 (a f F%diJ-N^lM{fi) 
sup {la JTnI ■ \\J\\l2„,, ■ ^16CIn(m) - N^In(m) 
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By the simple inequality for p,q > 0, 2pq — < p^, the last formula is bounded 
by a^\TN\ ■ II-/|I^2(t<;) ' i^^CIN^). We put a^N. On the other hand, 

lF~r^^y]v = ~ y, Vu,J{_oaixQ,a-iXQ)) ■ r]oa{xa.<T,xa) + on, 
ir^l ' - - 

where V„; is the directional derivative along m, and m, - ^ R''- 
By the entropy inequality, 

A^i^i;><logE;;exp|A^ + i/(P^|P;«). 

Since i/(PAi|P^) < |VV|C for some constant C, by the Feynman-Kac formula, we 
obtain ^ 

lim sup -^E^ r NF%dt < 16CT\\J\\l, , + \Vo\C' . 

N^oo \i-N\ Jo ' VI . J ; 

For a limit point of {Qn}n>\, Q*, we get 

Denote a countable dense subset of C''''([0, T] x T'') by J', we also get the following 
estimate: 



sup 



fl' 

Jo Jt' 



{-V u,J)pdndt - 16CT 



17 

Jo Jt" 



J djidt 



< \Vo\C'. 



See fT\ pp.107. Section 5.7 for details. Therefore for almost all p, there exists B{p) 
such that for every J e C°''([0, T] x T''), 

r r (-'^u,J)pdtidt - 16CT [ [ J^dpdt < B(p), 
Jo Jt' Jo Jt' 



that is, 



r r (-VuJ)pdiJdt <2Jl6CT f J^dpdt- ^[B{p). 

Jo Jt' ' \ Jt' 



This implies the linear- functional Ip : C"'\[0,T] x T'') ^ R defined by lp(J) := 
J^J-'VuJ)pdjj.dt is extended on L^([0, T] x T''). By the Riesz representation theo- 
rem, there exists V„_p € L?([0, T] x T'') such that 

r r {-"^ uJ)pdndt ^ r r jvupdndt 

Jo Jt'' Jo Jt' 

for every J e C''''([0, T] x T'') and every / = I, . . .,d. This yields Lemm dni □ 
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7.2 The uniqueness result 

We state the uniqueness result used in Sectio H?^ The following lemma follows from 
the argument by using the Grownwall inequality. 

Lemma 7.2. For any H € C''^([0, T] x T'), a weak solution of the quasi-linear partial 
differential equation 



is unique. 
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